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Abstract 

We give an exact solution to the generalized Langevin equation of motion of a charged Brow- 
nian particle in a uniform magnetic field that is driven internally by an exponentially-correlated 
stochastic force. A strong dissipation regime is described in which the ensemble-averaged fluctua- 
tions of the velocity exhibit transient oscillations that arise from memory effects. Also, we calculate 
generalized diffusion coefficients describing the transport of these particles and briefly discuss how 
they are affected by the magnetic field strength and correlation time. Our asymptotic results are 
extended to the general case of internal driving by correlated Gaussian stochastic forces with finite 
autocorrelation times. 
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The study of the Brownian motion of a charged particle in uniform and magnetostatic 
fields is of great importance in the description of the diffusion and transport of plasmas and 
heavy ions . In these references the respective authors modeled the 

dynamics of the charged particle through the use of a Langevin equation 1^ of the form 



Y{t)+jVit)-uV{t)x^ = riit) 



(1) 




where 7 is generally a constant dissipation rate tensor of rank two, u is the cyclotron 
frequency, B is the uniform magnetostatic induction field, and 77 (t) is the stochastic driving 
force per unit mass. Among the main objectives of these works are (i) to obtain expressions 
for the transport diffusion coefficient from the Langevin equation of motion in the presence 
„tdeUa-co.eIated (o. „M.e, and co^elated ,o. colo.d) „oi. 

(ii) the construction and solution of the associated Fokker-Planck equation [g], [?, 
drawback of these Langevin models emerges when the charged particle is driven by colored 
noise so that it generally does not achieve thermal equilibrium with the surrounding bath 
unless a restrictive condition is imposed on the dissipation rate and the autocorrelation time 
This situation is undesirable when the particle is driven internally (such as when the 
driving fluctuations are mainly due to collisions with particles constituting the fluid bath 
|lo|) and physical considerations require us to independently specify the dissipation rate 
and autocorrelation time. 

In this paper we address this problem by modeling the particle dynamics with a general- 
ized Langevin equation, replacing the dissipation force term of the Langevin equation with 
one that contains a memory integral and subsequently imposing a generalized fluctuation 



dissipation relation 12|, ll3|. In the presence of an exponentially-correlated stochastic driving 



force, this equation is solved exactly and representative ensemble-averaged trajectories given 
general initial conditions are presented. The average fluctuations in the position and velocity 
components are also calculated, allowing us to describe the thermalization and transport of 
these charged particles in uniform magnetic fields with the appropriate Gaussian noise and 
continuum approximations. Finally, we use Tauberian theorems to extend our asymptotic 
results for the case of arbitrary colored Gaussian driving noise with a finite correlation time. 

We consider a particle with charge q and mass m under the influence of a uniform magnetic 
induction field B and driven internally by a stochastic force per unit mass rjit). Since the 
Lorentz force on the particle is perpendicular to the induction field, the motion of the particle 
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in the direction of this field is governed by a free particle generalized Langevin equation with 
a well-known solution [l^. We therefore restrict the following analysis to the dynamics of 
the particle projected onto a plane perpendicular to the magnetic induction field B. Thus, 
the planar equation of motion expressed in Cartesian coordinates is 

v(t)+ /" r{t-t')v{t')dt' -iuj(TyV{t) = 'n{t), (2) 

Jo 

where r{t) is the memory dissipation kernel, u is the cyclotron frequency qB/m, and (Ty is 
the Pauli matrix ( ^ "^^ ) • 

We assume that the components of the stochastic force vector r]i{t) have Gaussian dis- 
tributions, vanishing average values {f]i{t)) = 0, vanishing cross-correlation functions, and 
identical exponential autocorrelation functions 

Here, 7 is the characteristic dissipation decay rate, r is the average correlation time, ks 
is Boltzmann's constant, and T is the absolute temperature of the bath surrounding the 
charged particle. These properties are based on the assumptions that (i) the bath is rota- 
tionally invariant about the field axis, (ii) particle collisions are accurately modeled as an 
Ornstein-Uhlenbeck process 

than the correlation time, a Fokker-Planck equation with a generalized diffusion coefficient 
proportional to the fiuctuations in position can be effectively constructed. 

To completely specify our problem, we finally apply the condition of internal driving by 
imposing a generalized fiuctuation-dissipation relation [i2\ on the dissipation memory kernel 
and find that 

T{t) = -e-*/^ t > 0. (4) 
r 

We begin our solution by obtaining the Laplace transform of the velocity 
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,^ [Vo + r7(a)], (5) 
s + V\ 

where I is an identity matrix, Vq is the initial velocity vector, and r(s) and f)(s) are the 
transforms of the dissipation memory kernel r(t) and the driving stochastic term rjit), 
respectively. The inverse matrix appearing in this transform may be formally expressed as 
a geometric series giving us 

V = — — f [^0 + ^]. (6) 
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FIG. 1: Poles of the transfer function described in the text for non-negative lot. Black dots represent 
the locations of these poles when wr = and arrowheads denote the direction of increasing lot. 



We split the resulting series into partial sums involving even and odd powers of <jy and use 

)btain 

1 + Ts) [(rs^ + s + 7)1 + iuj{\ + Ts)(T^ 



the fact that cr^ = I to obtain 



V 



[Vo + r7]. 



(rs2 + s + 7)2 + tu2(l + rs)2 
We identify the four poles of the transfer function multiplying [Vq + f\\ as 

— (1 + zu;r) ± a/1 — 47r — uj'^t'^ + i2a;r 



2t 



(7) 
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and their complex conjugates s*^. The locations of these poles in the complex plane are 
shown in Fig. [TJ A brief examination of this pole structure suggests that there exists a 
regime 7r < 1/4 where transient oscillations in the velocity components are negligible on 
the time scale of r for small cyclotron frequencies. Indeed, requiring the dominant relaxation 
time to be much smaller than the characteristic period of these oscillations under such a 
condition leads to the inequality <^ 7. On the other hand, when the product 7r is much 
larger than 1/4 these transient oscillations remain observable even as the cyclotron frequency 
vanishes, implying that they result from the memory of the system. This dramatic behavior 
will be demonstrated below as we construct the exact solutions to our generalized Langevin 
equation. 

With the poles of the transfer function identified, we can easily calculate the Laplace 
inverse of V by an integration over the appropriate Bromwich contour. The desired solution 
is 

V(t) = [(7i(t)I + ^^2(t)tT,]*[Vo5(t) + r7(t)], (9) 




FIG. 2: (Color online) Ensemble-averaged trajectories of internally driven charged particles in a 
uniform magnetic field directed in the z-direction. The product 7r is equal to 0.1 (a) and 1.0 (b). 
Adjacent symbols on each curve are separated in time by an interval r. 



where * is the Laplace convolution operator [17|, 6{t) is the usual Dirac delta function, and 
gi{t) and g2{t) are the real and imaginary parts, respectively, of the complex- valued function 



Bt l+iujT . ^ /3t 

cosh 1 smh — 

2r (3 2t 



(10) 



In the preceding expression, we have made the substitution (3 = a/1 — 47r — oo^t'^ + i2ujT. 
Since the stochastic force has a vanishing mean, we find that the ensemble-averaged velocity 
of the particle is 

(V(t))= [g,{t)l + g,{t)tcTy\Y,. (11) 

The average velocity asymptotically decays to zero as expected because of the presence of 
friction and the fact that the magnetic field does no work on the particle. The relaxation to 
equilibrium is generally damped biexponentially with characteristic times 2r/(l ± !He [/?]). 

By setting the origin of the Cartesian axes at the initial position of the particle and 
integrating (y{t)) we are able to calculate the average position of the Brownian particle as 

(X(t)) = [Gi{t)l + G2(t)2o-^] Vo, (12) 

where Gi{t) and G2{t) are the real and imaginary parts, respectively, of the function 

^(--^W- - cosh ^ - i^^I^sinh^l 



G{t) = 



e 2^ ^ 2^ 



(13) 
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Representative ensemble-averaged particle trajectories that were calculated from Eq. 
with the initial velocity Vq in the x-direction are depicted in Fig. [2l Consistent with the 
preceding analysis of the pole locations of the transfer function, we see that the average 



velocity components depict marked transient oscillations even in vanishing magnetic fields 
when the dissipation rate is high [ujt = 0.05 in Fig. [2]d). The long-time equilibrium position 
of the charged particle relative to its original position is 



7 ^ \ Vq 



= I ' I (14) 



As this expression is independent of the autocorrelation time, this result may be used to 
experimentally deduce the average dissipation rate 7 of appropriately prepared charged 
particle systems. 

The fluctuations in the position and velocity components may be expressed in terms of the 
functions g(t) and G(t) by noting the independence and equivalence of the noise components 



(Eq. [3]) and extending the approach outlined in Refs. [ij, [l5| yielding 



{{Av,r){t) = ^{l-\g{t)\'), (15) 



((AxO^)(t) 



m 



(^2^W[G(t')] dt' - \G{t)f^. (16) 



m 

Evaluating the remaining integral in the last equation leads to a lengthy expression, but the 
end result is a function consisting of a term linear in time plus transient terms, ((Axj)^) = 
2kBT'-ft /■m{uj'^ + '~f'^) + transient terms. This asymptotic linear dependence of the fluctuation 
in position on time tells us that the associated diffusion process is normal, as we discuss 
below. 

The fluctuation in velocity may be calculated from the mean-square velocity components 
(Eq. [T5|) and is graphed in Fig. [31 For the special case of a particle initially prepared with 
zero velocity, this quantity also corresponds to the average kinetic energy of the charged 



particle due to its projected motion on a plane perpendicular to the magnetic field |18 |. 
Since we have imposed the condition of internal driving and the mean velocity vanishes at 
long times, the quantity |m((AV)^) eventually relaxes to the expected equipartition energy 
ksT. In general, the typical time scale 6 of this relaxation is 

T 



which scales as 6' ~ co'^r^/7 for large cyclotron frequencies and small 7. Also, as we have 
discussed earlier, the approach of the velocity fluctuation to its equilibrium value for ■jt < 
1/4 and small values of ujt is essentially monotonic, while transient oscillations may be 
observed for small values of ut in the presence of strong friction (Figs. [3]d, [3}i). 




FIG. 3: (Color online) Time evolution of the velocity fluctuation ((AV)^) of the charged particle 
in units of its asymptotic value 2kBT/m. The product 7r is equal to 0.2 (a, c) and 5.0 (b, d). 

The effective one-dimensional transport diffusion coefficient Di{t) can be obtained from 
the ffuctuations in the components of the position vector (Eq. [T6!) through the definition 
Di(t) = J^((Axi)^). Calculating this derivative leads to 

2k T 

D.{t) = ^9\t{Gm~-g{t)]}, (18) 

where the overbar refers to the complex conjugate. In the limit of long times we find that 
these diffusion coefficients approach the constant value 



A~ = (19) 



which is independent of the correlation time r and correctly gives the zero-field value when 
a; — i> 0. This general result is reproduced by the analogous Langevin model {s], Q] only when 



the special condition 7r = 1 is imposed. Graphs displaying the transient characteristics 
of the diffusion coefficient are shown in Fig. HI In general, an increasing dissipation rate 
decreases the typical equilibration time as expected (Fig. Ht, Hji), although for large ujt 
the coefficient oscillates strongly before approaching its asymptotic value. The transient 
memory-induced oscillations of the diffusion coefficient in the strong friction regime are 
clearly seen when 7r = 5.0 (Fig. |Dd). 

We now turn our attention to generalizing our results to the case of a driving force with 
a finite characteristic correlation time but with an otherwise arbitrary autocorrelation func- 
tion. After imposing a fluctuation-dissipation condition as we have done above, Tauberian 
theorems allow us to extract the long-time behavior of such an autocorrelation function and 
the associated dissipation memory kernel by a leading-order truncation of the corresponding 
Laplace transforms [l6|. This approximation, which is valid for small s and large times, re- 
produces the same Laplace transform solution (Eq. [7j) upon a trivial rescaling of time. Thus, 
all the equilibrium values calculated here are generally applicable so long as the stochastic 
driving process is characterized by a finite autocorrelation time. 

In this paper we have analyzed the Brownian motion of a charged particle in the presence 
of a static and uniform magnetic field and driven internally by an exponentially-correlated 
stochastic force using a generalized Langevin model. The dynamics of the particle was de- 
scribed exactly when the driving noise is exponentially correlated in time and the velocity 
fluctuation and diffusion coefficient were calculated for all times. Of particular interest was 
a strong dissipation regime 7r ^ 1/4 that displayed damped oscillations in these quanti- 
ties even in the limit of vanishing cyclotron frequencies, which suggests that this transient 
behavior is memory- induced. Also, equilibrium values of these dynamical quantities were 
presented and demonstrated to be independent of the characteristic correlation time of the 
driving noise as required. With the use of Tauberian theorems we have shown that these 
equilibrium results do not depend on the specific form of the correlation function of the 
driving force so long as it possesses a finite autocorrelation time. 

J. P. E. acknowledges the support of a National Institute of Physics Faculty Grant. Part 
of this work was carried out by F. N. C. P. at The Abdus Salam International Centre for 
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FIG. 4: (Color online) Time evolution of the one-dimensional diffusion coefficient Di in units of its 
asymptotic value Df^ . The product 7r is equal to 0.2 (a, c) and 5.0 (b, d). 
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